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Abstract: Recently, the phenomenology of f\R) gravity has been scrutinized motivated 
by the possibility to account for the self-accelerated cosmic expansion without invoking 
dark energy sources. Besides, this kind of modified gravity is capable of addressing the 
dynamics of several self-gravitating systems alternatively to the presence of dark matter. It 
has been established that both metric and Palatini versions of these theories have interesting 
features but also manifest severe and different downsides. A hybrid combination of theories, 
containing elements from both these two formalisms, turns out to be also very successful 
accounting for the observed phenomenology and is able to avoid some drawbacks of the 
original approaches. This article reviews the formulation of this hybrid metric-Palatini 
approach and its main achievements in passing the local tests and in applications to 
astrophysical and cosmological scenarios, where it provides a unified approach to the 
problems of dark energy and dark matter. 
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1. Introduction 

One hundred years ago, Albert Einstein eompleted the mathematieal formulation of his revolutionary 
view of the gravitational interaetions in terms of eurved spaee-time. The spirit, eleganee, and 
experimental sueeesses of the original theory has eaptivated the international seientilie eommunity 
and the theory has been aeeepted as the standard model for gravity [1]. At its eentennial, pushed by 
new observational evidenees, the theory is at a dramatie erossroad. Its eontinuation as the referenee 
gravitational framework will imply that the universe is mainly eomposed by exotie sourees of matter 
and energy whose existenee is purely inferred from their gravitational effeets at the largest astrophysieal 
and eosmologieal seales [2-8]. However, if sueh sourees are not deteeted in any direet way, then we 
might be faeing a failure of one of the most original theories of the twentieth eentury. The important 
implieations of the two opposed alternatives, i.e., the seareh for unknown dark side eonstituents or the 
revision of gravitational theory, demand a eareful serutiny of the different possible seenarios. In this 
work, we eonsider the seeond ease, namely, the situation in whieh the gravitational dynamies may depart 
from that predieted by Einstein’s theory of General Relativity (GR) at ultraviolet and infrared seales. We 
mainly foeus on the latter regime. 

Given the sueeess of GR at relatively short seales (sueh as the Solar System, stellar models, or eompaet 
binary systems), the idea that modified dynamies eould arise at larger seales has been investigated in 
mueh detail over the last years. Theories in whieh the gravitational aetion eonsists of more general 
eombinations of eurvature invariants than the pure Einstein-Hilbert term have been investigated with 
speeial emphasis [9-17]. Erom these investigations it was soon notieed that the usual metrie formulation 
of alternative theories of gravity is generieally different from its Palatini (or metrie-affine) eounterpart 
(see [18] for a reeent review on the Palatini approaeh). Whereas the metrie approaeh typieally leads 
to higher-order derivative equations, in the Palatini formulation the resulting field equations are always 
seeond-order. The appealing eharaeter of the seeond-order equations of the Palatini formalism, however, 
is aeeompanied by eertain algebraie relations between the matter fields and the affine eonneetion, whieh 
is now determined by a set of equations eoupled to the matter fields and the metrie. The ease of f{R) 
theories is partieularly useful to illustrate the differenees between these two approaehes. In the metrie 
formulation, the objeet 0 = d//di? behaves as a dynamieal sealar field, whieh satisfies a seeond-order 
equation with self-interaetions that depend on the form of the Eagrangian f{R). In order to have an 
impaet at large astrophysieal and eosmologieal seales, the sealar field 0 should have a very low mass, 
implying a long interaetion range. It is well known, however, that light sealars do have an impaet at 
shorter seales, where their presenee is strongly eonstrained by laboratory and Solar System observations 
unless some kind of sereening meehanism is invoked [19-23]. In the Palatini ease, a sealar-tensor 
representation is also possible, but with the sealar field satisfying an algebraie rather than a differential 
equation. One then finds that the sealar 0 turns out to be an algebraie funetion of the traee of the 
stress-energy tensor of the matter, 0 = 0(T), whieh may lead, in models of late-time eosmie speed-up. 
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to undesired gradient instabilities at various eontexts, as has been shown by studies of eosmologieal 
perturbations [24,25] and atomie physios [26,27]. 

In this artiole we will review the hybrid variation of these theories, in whioh the (purely metrio) 
Einstein-Hilbert aetion is supplemented with (metrio-affine) oorreetion terms eonstruoted a la Palatini 
[28,29]. Given that metrie and Palatini f{R) theories allow the oonstruetion of simple extensions of GR 
with interesting properties and, at the same time, suffer from different types of drawbaeks, we initiated 
a program to establish bridges between these two seemingly disparate approaehes hoping to find ways 
to oure or improve their individual defioieneies. For that purpose, in a number of works we oonsidered 
a hybrid oombination of metrio and Palatini elements to oonstruot the gravity Lagrangian and found 
that viable models sharing properties of both formalisms are possible. An interesting aspeot of these 
theories is the possibility to generate long-range forees without entering into eonfliot with loeal tests of 
gravity and without invoking any kind of soreening meehanism (whioh would however require that at the 
present time the eosmologieal evolution reduees to GR.). The possibility of expressing these hybrid f{R) 
metrio-Palatini theories using a soalar-tensor representation simplifies the analysis of the field equations 
and the oonstruetion of solutions. In some sense, oonsidering a theory like R + f(jVj means that we 
retain all the positive results of GR, represented by the Einstein-Hilbert part of the aetion R, while the 
further “gravitational budget” is endowed in the metrio-affine f(JZ) eomponent. In faet it is well known 
that metrio-affine and purely metrie formalisms eoineide in GR, i.e., oonsidering the aetion R. On the 
oontrary, the two formalisms lead to different results oonsidering more generio funotions /{R) [18]. 

A related approaoh to study f{R) theories that interpolate between the metrio and Palatini families 
is that of the so-oalled C-theories proposed in [30,31]. There the spaeetime oonneotion is assooiated to 
the metrie = C{TZ)g^^ that is oonformally related to the spaeetime metrie but the relation may 
depend upon the a sealar eurvature TZ. This framework eontains the metrie, C{1Z) = 1, and the Palatini, 
C{TZ) = f'{R), formalisms as speeial limits, and one also finds that even when /{TZ) = IZ physieally 
distinot theories are possible. For further studies on variations of variational prinoiples see e.g. [32-38]. 

Other extensions of the f{R) framework modify the ooupling of matter to gravity by defining an aetion 
whioh depends linearly [39] or nonlinearly upon the matter Eagrangian [40-47], or its traee [48-50]. 
The new oouplings generally induee non-geodesie motion mediated by an extra foree orthogonal to the 
four-veloeity [51], whieh may have nontrivial effeets already in flat Minkowski spaee. Instabilities due 
to new nonlinear interaotions within the matter seotor are thus oommon in these theories [52,53]. We 
note, from this perspeotive, that in the hybrid metrio-Palatini approaoh, oonsidered in this artiole, it oan 
be expeeted that suoh instabilities in the matter seotor are absent, beeause the usual oonservation laws 
are satisfied. 

In this paper, we review the formulation and the main applieations of hybrid gravity models in 
late-time eosmologieal and astrophysioal soenarios. The artiole is organized in three main parts 
oonsidering the general formalism, the oosmology and the astrophysioal applieations. In Seotion 2, we 
start the disoussion oonsidering the aetion and the field equations of the hybrid metrie-Palatini formalism. 
In partieular, we disouss the soalar-tensor representation, the Cauehy problem, and more general hybrid 
theories than R + f{TZ). Seotion 3 is devoted to hybrid-gravity eosmology. We derive the Friedmann 
equations, eonstruet the related dynamieal system, and briefly eonsider some solutions. Furthermore, 
we analyse the eosmologieal perturbations in order to understand struoture formation in these theories. 
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We focus on the evolution of perturbations in the matter dominated era and vacuum fluctuations relevant 
to inflation. The weak field behaviour that is crucial for the Solar system precision tests of gravity 
is considered in Section 4, where we also discuss the galactic phenomenology and the astrophysical 
applications of hybrid gravity as an alternative to dark matter. In particular, we study the stellar dynamics 
and the theory of orbits, the generalization of virial theorem, the flat rotation curves of spiral galaxies, 
and the galactic clusters starting from the relativistic Boltzmann equation for collisionless systems of 
particles. The conclusions are drawn in Section 5. 

2. Hybrid metric-Palatini gravity: The general formalism 

Let us start our considerations by giving the basic features of the theory. In particular, we discuss the 
action and the field equations both in the so called f{X) and its equivalent scalar-tensor representations 
in both the Jordan and the Einstein frames. Next, the well-formulation and well-posedness of the Cauchy 
problem is considered. Finally, arbitrary hybrid gravity theories constructed from the metric and an 
independent connection are explored and it is shown that the special f{X) models avoid otherwise 
generic pathologies. 


2.1. Action and gravitational field equations 


The action is specified as [28,54] 

^ ^ [R + fin)] + , (1) 

where Sm is the matter action, = SttG, R is the Einstein-Hilbert term, R = g^''n^y is the Palatini 
curvature, defined in terms of an independent connection as 

K = = r - f+ f- f“/y), (d 

that generates the Ricci curvature tensor 7^^,^ as 


n = E" 


r“ +T 

fia,u ' 


■pA 

aA-*- {lu 


■pQ; pA 
^ au 


(3) 


Varying the action given by Eq. (1) with respect to the metric, one obtains the following gravitational 
field equation 

+ Fin)n^y — -fin)gf^i, = , (4) 

where the matter stress-energy tensor is defined as usual. 


, _ 2 di^Cm) 


(5) 


Varying the action with respect to the independent connection one then finds as the solution to the 
resulting equation of motion that is compatible with the metric Fin)g^u, conformally related to the 
physical metric g^y, with the conformal factor given by Fin) = dfin)/dn. This implies that 


n 


2^2(7^) 


Fin)^^^^'^'^'" 2 Fin) 


Fin),,Fin), y 


. ( 6 ) 
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The Palatini curvature, TZ, can be obtained from the trace of the field equation (4), which yields 

F{n)n-2f{n) = k^t + R = X. (V) 


Note that we can express IZ algebraically in terms of X if the form of f(JZ) allows analytic solutions. The 
variable X measures how much the theory deviates from the general relativity trace equation R = —n?T. 
These two observations shed light upon the structure of the theory and the central part that the ’’trace 
deviation” X plays in it. It is for this reason the R + f{TZ) hybrid metric-Palatini theories have also been 
called simply the ”f{X) theories”. 

Indeed, we can express the field equation (4) solely in terms of the metric and X as 


- F{X)R,, + ^ [F'(X)V„V“X + F"(X) [dXf] 




2 F{X) 


X^^X^y + K . 


( 8 ) 


Note that {dXY = X The trace of the field equations is now 


F'(X)V„V“X + 


F"(X) 


1 {F'{X)f 

2 F(X) 


{dXf + ^ [X + 2/(X) - F{X)R] = 0 , 


(9) 


while the relation between the metric scalar curvature R and the Palatini scalar curvature TZ is 


R(A') = fi+- 


'[F'ixy 

,V„V“F(X)' 

Knx), 

1 “ F(X) 


( 10 ) 


which can be obtained by contracting Eq. (6). 


2.2. Scalar-tensor representation 

Like in the pure metric and Palatini cases [55,56], the action (1) for the hybrid metric-Palatini theory 


can be turned into that of a scalar-tensor theory by introducing an auxiliary field A such that 

S = d^xV^[nAR + f{A) + fA{n - A)] + Sm , (11) 

where /a = df/dA and we have included a coupling constant for generality. Note that = 1 for 
the original hybrid metric-Palatini theory [28]. Rearranging the terms and defining (j) = /a, V{(p) = 
A/a — f{A), Eq. (11) becomes 

^ ^ [QaR + - 1/(0)] + Sm . (12) 

It is easy to see that this action is equivalent to our original starting point (1). Variation of the above 
action with respect to the metric, the scalar 0 and the connection leads to the field equations 

^AR^iu + (f^TZiMu — 2 d-~ 9fj.u = (13) 

TZ-V^ = 0, (14) 

= 0, (15) 
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respectively. 

The solution of Eq. (15) implies that the independent connection is the Levi-Civita connection of 
a metric This means that the relation (6) between the tensors and R^i, can be now 

rewritten as 

R-i.u = R^lu + ^df,(j)d,y(j) - ^ , (16) 

which can be used in the action (12) to get rid of the independent connection and to obtain the following 
scalar-tensor representation that belongs to the “Algebraic Family of Scalar-Tensor Theories” [57], so 
that one finally arrives at the following action 


5 




(f^A + (j))R + - V (0) 


+ S„ 


(17) 


It is interesting to point out that, by the substitution </> —>■ —(fi;0)^/6, the action (17) reduces to the 
well-known case of a conformally coupled scalar field with a self-interaction potential. Namely, this 
redefinition makes the kinetic term in the action (17) the standard one, and the action itself becomes that 
of a massive scalar-field conformally coupled to the Einstein gravity. Of course, it is not the Brans-Dicke 
gravity where the scalar field is massless but now we have a nonzero V(</>) as given in (12). 

In the limit 0, the theory (17) presents the Palatini-/(7?.) gravity, and in the limit —)■ oo 

the metric f(R) gravity [57]. Apart from these singular cases, the more generic theories with a finite Oyi 
thus lie in the ’’hybrid” regime, which from this perspective provides a unique interpolation between the 
two a priori completely distinct classes of gravity theories. In fact, we have arrived at Brans-Dicke type 
of theories specified by the non-trivial coupling function 


30 

“ 20 - 20a ’ 


(18) 


that generalises the ujbd = 0 and cjbd = —3/2 cases, which correspond to the scalar-tensor 
representations of the metric f(R) and the Palatini-/(7^) gravities [14], respectively. 

Using Eq. (16) and Eq. (14) in Eq. (13), the metric field equation can be written as 


(Ua + 4’)R^lu 




+ 2^^*^ Vq,V“ 0) -f 'VijR7 


3 


9 ^ 09^0 , 


(19) 


or equivalently as 

+ 0)67^1/ = + 'Vu4> — — ^V^0Vi/0 -f —Va0V'^05'^j/ — (20) 

from which it is seen that the spacetime curvature is generated by both the matter and the scalar field. 
The scalar field equation can be manipulated in two different ways that illustrate further how the hybrid 
models combine physical features of the ojbd = 0 and ojbd = —3/2 scalar-tensor models. 

First, tracing Eq. (13) with g^'^, we find —— (pTZ -\- 2V = k^T, and using Eq. (14), it takes the 
following form: 

2V - 014 = k^T + QaR ■ (21) 

Similarly as in the Palatini (ujbd = —3/2) case, this equation tells us that the field 0 can be expressed as 
an algebraic function of the scalar X = n?T -f VLaR, i.e., 0 = 0(A). In the pure Palatini case, however. 
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0 is just a function of T. The right-hand side of Eq. (19), therefore, besides eontaining new matter terms 
assoeiated with the traee T and its derivatives, also eontains the eurvature R and its derivatives. Thus, 
this theory ean be seen as a higher-derivative theory in both matter and metrie fields. However, sueh an 
interpretation ean be avoided if R is replaeed in Eq. (21) with the relation 

R = n + I (22) 


together with TZ = V^. One then finds that the sealar field is governed by the seeond-order evolution 
equation that beeomes, when = 1, 




(1 -f 0)V0] _ 0K" 


T, 


(23) 


whieh is an effeetive Klein-Gordon equation. This last expression shows that, unlike in the Palatini 
{ojbd = —3/2) ease, the sealar field is dynamieal. The theory is therefore not affeeted by the mieroseopie 
instabilities that arise in Palatini models with infrared eorreetions [18]. 

Einally, we ean make a eonformal transformation into the Einstein frame of these theories. The 
eonformal resealing we need to aehieve this is given by 


Qfiv — (0 T 1 


(24) 


and the Einstein frame Eagrangian then beeomes 

C = R I _ ^(0) 

20 (0 -f O^) (0 -f Oa) 


(25) 


This ean be further put into its eanonieal form by introdueing the resealed field 0 as 


0 = Oa tan 


2 



(26) 


The vaeuum theory then beeomes a eanonieal sealar theory with a very speeifie potential (stemming of 
eourse from the original funetion f{7l)) in the Einstein frame. 


2.3. The Cauchy problem 


The dynamieal equivalenee with sealar-tensor theories shown above is useful to diseuss the 
well-posedness of the Cauehy problem for hybrid /(X)-gravity in vaeuo and eoupled to standard matter 
sourees. Eor previous studies of the Cauehy problem in different formulations of f{R) theories, see 
[58-62]. 

We begin by proving the well-posedness of the Cauehy problem in vaeuo, making use of the 
equivalent formulation (20) and (23). As we shall argue, the same eonelusions hold in the presenee 
of standard matter sourees satisfying the usual eonservation laws V^T^jy = 0. 

Borrowing definitions and notations from [63], the key point of our diseussion is the introduetion of 
suitable generalized harmonie eoordinates, defined by the eonditions 




1 

(! + </>) 


V'^0. 


with := : 


(27) 
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As we shall see, the gauge (27) allows us to develop a second order analysis very similar to the one used 
in GR [63]. We notice that the generalized harmonic gauge (27) is a particular case of the one introduced 
in [64] to prove the well-posedness of the Cauchy problem for a certain class of scalar-tensor theories of 
gravity. 

Let us start with rewriting Eqs. (20) in the form [65] 


R 




(! + '#>) L 





(28) 


where 

■= , (29) 

plays the role of an effective energy-momentum tensor. We recall that the Ricci tensor can be expressed 
as [63] 

= RU + \ [%<rdu {f; + H-) + {f; + h^)]^ (30) 

with 

Rtu ■= - + Af,^{g, dg ), (31) 

where only first order derivatives appear in the functions Due to the assumed gauge condition 
= 0 and the explicit expression of H^, from (30) and (31), we get the following representation 


Ri.v = + B/iAg, 0 , dg, d(j)), (32) 

where the functions depend on the metric g, the scalar field 0 and their first order derivatives. At the 
same time, using Eq. (23) to replace all terms depending on the divergence V ^(p, the right hand 

side of (28) can be expressed as 


(1 + ■« i 




(! + </>) 


dlA + Cf,^{g,(j),dg,d(j)), 


(33) 


where, again, the functions 67^1, depend only on first order derivatives. A direct comparison of Eq. (32) 
with Eq. (33) shows that, in the considered gauge, Eq. (28) assumes the form 


= Df,^{g, (p, dg, dcp). 


(34) 


The conclusion follows that Eq. (23) together with Eq. (34), form a quasi-diagonal, quasi-linear 
second-order system of partial differential equations, for which well known theorems by Eeray [63,66, 
67] hold. Given initial data on a space-like surface, the associated Cauchy problem is then well-posed in 
suitable Sobolev spaces [63]. Of course, the initial data have to satisfy the gauge conditions = 0 as 
well as the Hamiltonian and momentum constraints 

/^ = 0,...,3, (35) 

{l + (p) 

on the initial space-like surface. In connection with this, we notice that, from Eq. (23), we can derive 
the expression of the second partial derivative d^cp and replace it on the right hand side of (35), and 
thus obtaining constraints involving no higher than first order partial derivatives with respect to the 
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time variable To eonelude, we have to prove that the gauge eonditions = 0 are preserved in a 
neighbourhood of the initial spaee-like surfaee. To this end, we first verify that the divergenee of the 
gravitational field equation (20) vanishes, namely 

Vq(l + 0)G'^.-S^,] = 0. (36) 

Taking into account the identities = 0 and (V^(j)) — Vi^V^V^) 0, 

automatically satisfied by the Einstein and Ricci tensors, we have 

[(1 + 0)G^. - + \v{cl>)g,)j . (37) 

On the other hand, inserting the content of Eq. (21) (in this case, with T = 0) into the trace of the field 
equation (20), we end up with the identity 

The identities (36) follow then from a direct comparison of (37) with (38). 

Now, if and 0 solve the reduced field Eq. (34) and the scalar field Eq. (23), then we have 

(1 + 4,)G'“ - S'*" = -T±i) (gi“d,F; + . (39) 

Identities (36) imply then that the functions satisfy necessarily a linear homogeneous system of wave 
equations of the form 

= 0 , (40) 

where F*^ are known functions on the space-time. Since the constraints (35) amount to the condition 
(OqF^ = 0 on the initial space-like surface, a well known uniqueness theorem for differential systems 
such as Eq. (40) assures that F^ = 0 in the region where solutions of Eqs. (23) and (34) exist (see also 
[63]). 

The illustrated analysis also applies in the case of couplings to standard matter sources such as 
electromagnetic or Yang-Mills fields, (charged) perfect fluid, (charged) dust, Klein-Gordon scalar fields 
[61], so showing the well-posedness of the Cauchy problem for /(X)-gravity in presence of standard 
matter fields. Indeed, when matter sources are present, Eqs. (23) and (34) have to be coupled with the 
matter field equations. Applying the same arguments developed for GR [63,68,69], it is easily seen that, 
in the generalized harmonic gauge (27), the matter field equations together with Eqs. (23) and (34) form 
a Eeray hyperbolic and a causal differential system admitting a well-posed Cauchy problem [67]. In 
addition to the well-known results by Bruhat, the crucial point is again that the field equations of matter 
field imply the standard conservation laws = 0 [39]. In summary, the hybrid metric-Palatini 

gravity satisfies the well-formulation and well-posedness of Cauchy problem for standard forms of matter 
and then, in this sense, it is a viable theory. 
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2.4. More general hybrid metric-Palatini theories 

The “hybrid” theory space is a priori large. In addition to the metric and its Levi-Civita connection, 
one also has an additional independent connection as a building block to construct curvature invariants 
from. Thus one can consider various new terms such as 

R^'^R^u, etc. (41) 

Though an exhaustive analysis of such hybrid theories has not been performed, there is some evidence 
that the so called f{X) class of theories we are focusing our attention upon here is a unique class of viable 
higher order hybrid gravity theories. In the more restricted framework of purely metric theories, it is 
well known that the f{R) class of theories is exceptional by avoiding the otherwise generic Ostrogradski 
instabilities by allowing a separation of the additional degrees of freedom into a harmless scalar degree of 
freedom [70]: as we have already seen, such a separation is possible also for our hybrid-/(X) theories. 
Furthermore, it turns out that this feature is a similar exception in the larger space of metric-affine 
theories, since a generic theory there is inhabited by ghosts, superluminalities or other unphysical degrees 
of freedom. 

As a representative class of more general theories, actions of the form 

^ ^ f{R, Qh) , Qh = R^'^K. (42) 

were studied by Tamanini [71], who determined the precise field content of this action in the weak-field 
limit. Variation of (42) with respect to the metric produces the following field equations: 

I^rRiiu — — Xfj,Vuf,R + f,'R.Riiu + ‘^f^qR^RuX 

+^n(/eV) + ^9,.v„v^(/gB“'=)-VAV,„(/4fl;,) = (43) 

where f^R, f^-ji and f q axe the derivatives of / with respect to R, TZ and Qh respectively. The solution to 
the equation of motion for the connection on the other hand dictates that it is the Levi-Civita connection 
of the metric given by 

^ , where + f qR^'^ ■ (44) 

V~9 

Using this one can eliminate the auxiliary metric in terms of the physical metric g^,^. 

Considering perturbations = g^^ — around Minkowski space g^y = g^y, and inverting the 
linearised field equations for the physical metric then gives us the propagators for the graviton and the 
additional degrees of freedom that may be present in h^j^y. The propagator is defined by 

= (45) 

where Taji represents the linearised stress energy source. In the formalism of the spin-projector operators 
employed in Ref. [72] and more pedagogically reviewed in Ref. [73], the result can be given in Fourier 
space (where basically □ —— A:^) in terms of two functions a and c as 

•p2 -pO 

l2tt _ __ 

a(_p) _3c(_A;2) 


( 46 ) 





Universe 2015 , xx 


11 


where pieks up the spin-2 and the sealar modes of the fluetuations, see Refs. [71-73] for 
details. The funetions a and c ean be determined immediately given a theory of the form (42). They 
depend upon the eombinations 


in the following way: 


A = 


“(□) = /S’+ 


6/SS+/:”’ 


2/ 


(0) 

,TZ 


and B = 


f 


( 0 ) 

,Q 


/: 


( 0 ) ’ 
Tl 


(47) 


' ,R 


,7^ 


,Q 4 


c(n) = /S’+/S’- 2 (/SI + 4/SS +/f) ° + 


/SS (6.4 + B) + / 


( 0 ) 

,Q 


, B 
2A+- 


(48) 


(49) 

(50) 


Let us then enumerate some speeial eases. To simplify things we assume = 0. 
2.4.1. Metrie f{R) models 

In the pure metrie f{R) ease, = ^4 = 0 and we have 


_ -rT-a/37<5 I 
“ ^'-GR “T 


2 R- + (3/SS)-‘ 


(51) 


Thus we have an extra sealar degree of freedom, as we expeet sinee the f{R) models are known to be 
equivalent to Brans-Dieke theories with a vanishing parameter ujbd = 0. The mass of the “sealaron” is 
and as long as f”{R) > 0 the theory is stable, otherwise a taehyonie mass spoils the 
stability around Minkowski spaee. 

2.4.2. Palatini f{TZ) models 

As already diseussed, the Palatini-type f(JZ) models are equivalent to Brans-Dieke theories with 
the parameter ujbd = —3/2. This partieular value eorresponds to vanishing kinetic term of the field, 
which thus is nondynamical. Therefore we expect that no additional scalar degree of freedom should 
appear. For a proper normalisation we may assume that = 1, and we have now of course that 

f%R = 

(52) 


_ tto/37(5 


confirming our expectation. 

2.4.3. Hybrid f{X) models 

It was already remarked in [57], in Ricci-flat spacetimes the f{X) theories share the properties of 
Palatini-/(77.) theories, which in vacuum reduce to GR with a possible cosmological constant. Therefore 
it is not a surprise that we find no new propagating degrees of freedom in Minkowski vacuum. 


_ tto / 37(5 
‘■‘■fix) ~ ‘■‘■GR 


( 53 ) 
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Interestingly though, this elass of theories is not equivalent to either of the previous two oases, sinee 
when one oonsiders eurved spaoetimes, a new soalar degree of freedom appears. In this sense, the f(X) 
gravity is a quite minimalistie soalar-tensor extension of GR, as the soalar propagates only in the presenoe 
of baokground eurvature. 

2.4.4. The hybrid f(R, 7Z) models 


The generalized hybrid Rioei soalar theories were introdueed in [74,75] and found to have 
qualitatively different properties eompared to the more restrioted elass of /(X) models desoribed above. 
In partioular, the f{R,TZ) were shown to be equivalent to a elass of bisoalar-tensor theories. These 
theories have an extra propagator with a double pole, oorresponding to two propagating 

sealar degrees of freedom. We ean easily deduoe the masses of these soalar fields. We get 



where we have defined for oonvenienoe 


S 



(54) 


(55) 


We note that the sealar partiele with mass squared m?_ oorresponds to the sealaron appearing in (51) in 
the limit of pure f{R) gravity, but in general now has a shifted mass. The other sealar is a new partiele 
that ooours due to nontrivial dependenoe upon TZ, and unlike in the ease of /(X) gravity, it propagates 
also in Riooi-flat spaoes. The oondition that neither of the soalars has a taohyonio instability, is given by 


fjl > 0 ’ f^RR + ^f^RTl S > 0 . 


(56) 


The residues at the two poles oorresponding to these masses are 


r± = 



(57) 


In order for neither of these sealars to be a ghost, we should have both r+ > 0 and r_ > 0. The seoond 
oondition would require that 

f^RR + — *5' < 0, (58) 

in oontradiotion with (56). It seems then that we oannot avoid both taohyons and ghosts in this theory. 


2.4.5. The hybrid Riooi-squared f{Tl, Q) theories 

Let us finally eonsider the Qjr-invariant. For simplieity, we restriet to models here without nonlinear 
dependenoe on the metrie Rieei soalar; it is easy to see that this does not affeet our eonolusions essentially. 
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Basically the graviton propagator acquires its structure from the function a(n) in (48), and now only the 
higher-derivative term Qh modifies it. We can arrange the result for the propagator in the form 



The sixth order theory we have at hand has a modulated graviton propagator which adds two extra poles. 
In addition, there appears a scalar propagator that has five poles. This is in a quite drastic contrast with 
respect to the metric Q-theory which contains only one additional spin-2 particle and features fourth 
order field equations. We need not analyze in detail the properties of the new degrees of freedom here, 
since it is obvious the theory as such is seriously haunted by ghosts and thus not physical. It is easy to 
convince oneself that this occurs very generically once one builds the action from any hybrid curvature 
invariant - with the exception of IZ in the specific case of separable functional dependence R + f{TZ)- 
These considerations corroborate our claim that the f{X) theories are of special theoretical interest. 
In the rest of the review, we shall discuss their phenomenology. 


3. Hybrid-gravity cosmology 


In order to study the cosmology of the metric-Palatini theories, we choose in this section to employ 
the scalar-tensor formulation derived above (17): 





(flyi -f (j))R -\- 


3 


{d<pf 


2kV(0) 


+ Sr, 


(60) 


where 

= \[r{(p)(p-f{r{(j)))], r{(j)) = . (61) 

Cosmology of the equivalent theories has been investigated also using the f{X) formulation in terms 
of the purely metric quantities, and we refer the reader to the recent phase space analysis for the most 
complete global analysis of the cosmological dynamics of these theories [76]. Here we will first write 
down the cosmological equations in the formulation (60) and then have a brief look at the phase space 
of exact solutions for these equations. Then we will analyse the formation of cosmological large-scale 
structure in these models. 


3.1. Background expansion 

The flat Friedmann-Robertson-Walker metric is defined as 

ds^ = —-f a^{t) (dx^ -f dy^ + , (62) 

where the rate of time-evolution of the scale factor a{t) is conveniently parameterised by the Hubble 
rate H{t) = {da{t)/dt)/a{t). In the following we will not write explicitly the time dependence of the 
cosmological background quantities, and denote time derivative by an overdot. This eases the notation 
and we can write for example H = d/a. In the following we will mainly be interested in accelerating 
dark energy -like dynamics; for a study of Einstein static spaces, see Ref. [77]. 
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3.1.1. The Friedmann equations 

The Friedmann equations that govern the evolution of H ean always be written in terms of the effeetive 
energy density and pressure, respeetively defined as 

= K^PeS, (63) 

^2 

H = -y (Peff+Peff) , (64) 

where for the theory defined by the aetion (60) we obtain the following effeetive souree terms 

(^A + 0) (0) - 3iT0 + K^Pm , (65) 

A(p 

{Ha + 0) K^PeS = --^0^ - (0) + 0 + 2i?0 + K^Pm , (66) 

A(p 

respeetively. The eonservation equations for the matter eomponent and the sealar field are 

pm + ^H{pra+Pm) = 0, (67) 

0 + 3770-^ +j0/2-^«W'(0) = 0. (68) 

2(p 6 6 

Reealling that R = 6(277^ + 77) and using Eqs. (63) and (64), we ean rewrite the Klein-Gordon equation 
as 


where for notational simplieity, 77'(0) is defined by 

2k^0 


{Pm ^Pm) 0 ; 


f/'(0) = yy [24^(0) - ((7^ + 0) 4^'(0)] 


(69) 


(70) 


As a eonsisteney eheek one ean verify that the Klein-Gordon equation together with the matter 
eonservation allows to derive (64) from (63). By eombining Eqs. (68) and (69), we find that 


21/(0) - 4"'(0)0 = ^ {HaR + K^Tm) = 


(71) 


The solution for 0 = 0(X = 0) gives us the natural initial eondition for the field in the early universe. 
The asymptotie value of the field in the far future may then be dedueed by studying the minima of the 
funetion U (0) defined by Eq. (70). 

3.1.2. Dynamieal system analysis 


Cosmologieal dynamies ean be addressed by taking into aeeount a suitable dynamieal system. Eet us 
introduee the dimensionless variables 

3^2 ’ ^ ’ y-x,N, ^ “ 3 ^ ’ 

where N = log a is the e-folding time. The Friedmann equation (63) can then be rewritten as 


— 




( 72 ) 






Universe 2015 , xx 


15 


Due to this constraint, the number of independent degrees of freedom is three instead of four. We choose 
to span our phase space by the triplet {x, y, z}. The autonomous system of equations for them reads as 


X,N 

= y, 






(74) 

y,N 

_ 2x + 2/ J 

[ (3Wm - 

-i)y^ 

' + Ax [{3Wm -1)2/- 

-3(1 + Wm) z] 



—4x^(1 

3'U^^ 

-2u{ 

x)z) + AVLa 

[3X {Wm 

-l)y + y^ -x^(2- 6wm - Au{x)z)] 

[ ,(V5) 


= " \ 
AVLax 1 

{^WrfTl 

i)y^ 

+ Ax [{3Wm 

-^)y- 

■3(1 + Wm) z] 



+AVLax 

(3 + ?>w, 

m,R U 

{x)y) + 4x^ 

{3Wm 

1 + 2u{x)z) 1 . 

(76) 


We have defined u{x) = V'{(f))/V{(j)), that becomes a constant iff the potential is exponential. The 
relevant fixed points appear in this system. In particular, we have the matter dominated fixed point 
where x = y = z = 0 and Wes = Wm, and the de Sitter fixed point' that is described by Wes = — 1 and 

x^ = {2— VLau^)/u* ■, r/* = 0, z^ = 2/u^. (77) 

We denote the asymptotic values corresponding to this fixed point by a subscript star. In particular, the 
asymptotic value of the field x* is solved from the first equation in (77) once the form of the potential is 
given. As expected, this value corresponds to minimum of the effective potential (70), U'{x^) = 0. To 
construct a viable model, the potential should be such that we meet the two requirements: 

• The matter dominated fixed point should be a saddle point, the de Sitter fixed point an attractor. 
Then we naturally obtain a transition to acceleration following standard cosmological evolution. 

• At the present epoch the field value should be sufficiently close to zero. Then we avoid conflict 
with the Solar System tests of gravity (this will be clarified in Section 4.1). 

Note that the simplest metric f{R) theories that provide acceleration fail in both predicting a viable 
structure formation era and the Solar system as we observe it. The Palatini-/(7?.) models on the other 
hand can be ruled out as dark energy alternative by considering their structure formation or implications 
to microphysics, if such a theory is regarded consistent in the first place. As shown here and explored 
further below, f{X) gravity models exist that are free of these problems. 

To summarize: the field goes from to /*, where the former is given by 21/(0*) = l/'(0j)0j 
and the latter by 2V (0*) = (D^ + 0*)l/'(0*)- We just need a suitable function V (0), i.e. f{Tl) in such 
a way that the slope will be downwards and 0* near the origin. 

We refer the reader to [76] for a more complete and detailed phase space analysis of the cosmological 
background dynamics. 


1 


In addition, there exists the fixed point x 


flA corresponding to some kind of singular evolution. 
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3.1.3. On cosmological solutions 


As a specific simple example, let us eonsider in more detail the speeifie ease of de Sitter solution in 
vaeuum when VLa = 1. Then the modified Friedmann equations take the form 


?>H‘^ 

2H 


1 

1 + 0 

1 

1 + 0 


- 3 ^ , 

3 0^ 

-K‘^{pm+Pm) + H^+ -0 


5 


and the sealar field equation (23) beeomes 

0 + 3Ff0 - ^ + ^[21^ - (1 + <P)V^] = -^{Pm - Spm). 


(78) 

(79) 


(80) 


To further speeify the set-up, eonsider a model that arises by demanding that matter and eurvature satisfy 
the same relation as in GR. Taking 

V^(0) = K) + 100^ (81) 

the traee equation automatieally implies R = —k^T + 2Vq [28,54]. As T —)■ 0 with the eosmie 
expansion, this model naturally evolves into a de Sitter phase, whieh requires Vq ~ A for eonsisteney 
with observations. If Vi is positive, the de Sitter regime represents the minimum of the potential. The 
effeetive mass for loeal experiments, = 2(Vo — 21/i0)/3, is then positive and small as long as 
0 < Vq/Vi. For suffieiently large Vi one ean make the field amplitude small enough to be in agreement 
with Solar System tests. It is interesting that the exaet de Sitter solution is eompatible with dynamies of 
the sealar field in this model. 

The aeeelerating dynamies that drive the f{X) theory towards its general relativistie limits today have 
indeed been realised in several speeifie models [28,54,78,79]. Our preliminary phase space analysis 
eonfirmed the existenee of de Sitter attraetor solutions, and the reeent study of eosmology in terms of 
dynamieal system analysis extends this result to more general models [76]. Analytic solutions were 
presented also in Ref. [54] as well as in Ref. [80], there using a Nother symmetry teehnique. A designer 
approaeh was devieed by Lima [78] to reeonstruet preeisely the standard ACDM expansion history by 
a nontrivial f{X) model, and finally, two families of models were eonstrained by eonfronting their 
predietions with a eombination of eosmie mierowave baekground, supemovae la and baryonie aeoustie 
oscillations background data [79]. 


3.2. Cosmological perturbations 

To understand the implications of these models to the eosmologieal strueture formation, we will derive 
the perturbation equations and analyse them in some speeifie cases of interest. This paves the way for a 
detailed eomparison of the predietions with the eosmologieal data on large seale strueture and the eosmie 
mierowave baekground. For generality, we will keep the parameter Ha in the formulas in this seetion. 









Universe 2015 , xx 


17 


3.2.1. Field equations and conservation laws 


We work in the Newtonian gauge [81], which can be parameterized by the two gravitational potentials 
$ and 'k, 

= - (1 + 2^) + a^(t) (1 + 2$) df^ . (82) 

As matter source we consider a perfect fluid, with the background equation of state w and with density 
perturbation 5 = 5pm/Pm, pressure perturbation 5pm = c^s^Pm and velocity perturbation v. 

The 0-0 part of the field equations is 


^$ + 3 IF- 




H(j) 


^ 2 (Oa J y + 0 40 + 0) 




2f2A + 0 

where we have denoted p = 50. The Raychaudhuri equation for the perturbations reads 


(83) 


Q[H'^ + 2H] -2^ + 


6 


Oa + 0 




^ + 3 2i7 - 




VLa + 0 


$ - 'k - 6<l> = 


20 


+ 0 

The 0-i equation is 

-in. 


^ i^^pm H“ + I QH + QH + 3— — H—- j + 3 [ -— j (p + 3^^ 


.(84) 


<k + <k = 


{pm + Pm) aVm + \ H + ^ ] p + If 


2 ((2a + 0) y 2 ((2 a + 0) 

Note that the set of perturbed field equations is completed by the off-diagonal spatial piece: 

\k -f <k = -- 


(85) 




(2a + 0 

Assuming a perfect fluid, the continuity and Euler equations for the matter component are 

S + 3H — w) S = — (1 -f w) ^36 — —V 


( 86 ) 


i; -I- (1 — 3c^) Hv = - I \k -f 


-5 


1 -f te 


(87) 

( 88 ) 


respectively. The linear part of the Klein-Gordon equation is then compatible with the above system. 
For completeness, it is 


This completes the presentation of the field equations and the conservation laws. For the equations in 
the synchronous gauge, see [78]. 
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3.2.2. Matter dominated eosmology 


Let us eonsider the formation of strueture in the matter-dominated universe, where te = = 0. 

In this subseetion, we shall eonsider seales deep inside the Hubble radius. This so ealled quasi-statie 
approximation is well known in the literature and indeed we will arrive at a similar result as have 
been known to apply for sealar-tensor theories sinee early studies [82] that have been more reeently 
generalised to wide a variety of eoupled dark seetor models [83]. The approximation negleets a 
fluetuating degree of freedom that is expeeted to be insignifieant at small enough seales. For analyses 
of the applieability of the approximation with different assumptions on the eosmologieal models, see 
[83-87]. 

In the quasi-statie subhorison limit the spatial gradients are more important than the time derivatives 
and, eonsequently, the matter density perturbations are mueh stronger than the gravitational potentials. 
Combining the eontinuity and the Euler equation in this approximation, one obtains 

6 = -2H6-—^. (89) 

We need then to solve the gravitational potential. Let us define 11 = a‘^pm5/k‘^ and write the field 
equations and the Klein-Gordon equation at this limit in a very simple way as 


+ = n-v?, (90) 

{Qa + <P){^ + ^) = (91) 

—2 (O^ -f 0) \I/ = n -|- (^ , (92) 

= -20 (^ + 2$). (93) 


We immediately see that one of the equations is (as expeeted) redundant, and that the fh is (as usual) 
proportional to IT, where now the proportionality is given as a funetion of the field 0. Our result is 

6 + 2H6 = 47rG,sPmS, (94) 

with 






-G. 


(95) 


Oa -f 0) 

This shows that instabilities ean be avoided in the evolution of the matter inhomogeneities, in eontrast to 
the Palatini-/(7^) models and some matter-eoupled sealar field models (reeall our theory ean be mapped 
into sueh in the Einstein frame). Equation (94) provides a very simple approximation to traek the growth 
of strueture aeeurately within the linear regime during matter dominated eosmology. Confrontations of 
speeifie model predietions with the present large seale strueture data and foreeasts for the eonstraints 
from future experiments, in partieular the Euelid mission [110], is interesting work to be done. 


3.2.3. Vaeuum fiuetuations 


The propagation of our sealar degree of freedom in vaeuum is also a erueial eonsisteney eheek on the 
theory. Let us set = 0. Let us eonsider the eurvature perturbation in the uniform-field gauge In 
terms of the Newtonian gauge perturbations this is 

C = $ - 

0 


(96) 
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After somewhat more tedious algebra than in the previous ease, we obtain the exact (linear) evolution 
equation 


c+ 


3ff-2 


4 > + ‘2H + 0) 


(j) + 2H + 0) 


0 (Oa + 4 >) ^ 200 0^ (flA + 0)^ 0 


02 l0(flyi + 0) 1 _ 03 ^ 


C = -^C. 

/-r^ 


(97) 

The friction term depends on the perturbation variable we consider, but the perturbations at small scales 
still propagate with the speed of light, as in canonical scalar field theory. This excludes also gradient 
and tachyon instabilities in the graviscalar sector. Now equation (97) can be used to study generation 
of fluctuations in /(X)-inflation. Construction of specific models and their observational tests are left 
for forthcoming studies; the Einstein-frame formulation (25) might present a convenient starting point 
for that as it, given the function f(jVj, presents directly the relevant inflationary potential in terms of the 
canonic field. 


4. Astrophysical applications 

In this Section we develop methods to study dark matter phenomenology in the hybrid models; for 
studies of dark matter generated by metric f{R) modification, see e.g. [88,89]. 

Hybrid gravity allows to address several issues related to dark matter dynamics ranging from galaxies 
to galaxy clusters. These self-gravitating structures can be probed by studying the motion of test particles 
(stars for galaxies and galaxies themselves for galaxy clusters) moving into a gravitational potential. 
The behaviour of rotational and dispersion velocities of such test particles can be explained within the 
framework of the gravitational potential derived from the theory. For example, the tangential velocity 
can be explicitly obtained as a function of the scalar field of the equivalent scalar-tensor description. 
The model predictions can be compared with samples of rotation curves of spiral galaxies and galaxy 
clusters, respectively. The possibility of constraining the form of the scalar field and the parameters 
of the model by using the stellar velocity dispersions is also analysed. Furthermore, the Doppler 
velocity shifts are also obtained in terms of the scalar field. Finally suitable generalizations of the virial 
theorem and the relativistic Boltzmann equation allow to construct a self-consistent theory for galaxy 
clusters. In conclusion, all the physical and geometrical quantities and the numerical parameters in the 
hybrid metric-Palatini model can be expressed in terms of observable/measurable quantities, such as the 
tangential velocity, the baryonic mass of the galaxy, the Doppler frequency shifts, the dispersion velocity, 
the geometrical quantities characterizing the clusters of galaxies respectively. These results open the 
possibility of testing the hybrid metric-Palatini gravitational models at the galactic or extra-galactic scale 
by using direct astronomical and astrophysical observations. 

Fet us start our considerations by dealing with the weak field limit of the theory. 

4.1. The weak field limit 

It is of paramount importance to determine the post-Newtonian parameters of the theory as they 
determine the compatibility of the theory with the local precision gravity tests. For post-Newtonian 
analysis of the metric and Palatini f{R) theories, see e.g. Refs. [55,56], and for a unified analysis Ref. 
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[90]. Here, in partieular, we are interested in the parameter 7 that is basieally the fraetional difference of 
the Newtonian potentials in (82) in the limit where the cosmological expansion can be neglected, a = 1. 

To this end, we need to consider the perturbations of Eqs. (19) and (23) in a Minkowskian background. 
The usual procedure is to assume cf) = (f)Q where 0o is the asymptotic value of the field far away 

from the local system (and should be given by the cosmological background solution), and to take a 
quasi-Minkowskian coordinate system in which rj^^, + h^y, with \h^y\ 1. Then it is easy to see 

that the quadratic terms and in Eqs. (19) and (23) do not contribute to the linear order. 

The potential terms in Eq. (23) can be linearized as follows (in this subsection we denote = E'(0)): 


t[2V - (1 + 0)E^] ^ ^[2Eo - (1 + 0o)^^o] - I [0o(l + - 2Vo] • (98) 


The zeroth-order term in this equation is due to the background, and can be absorbed into a coordinate 


redefinition. The coefficient of the first-order term can be interpreted as a mass squared. The linearized 
scalar field equation is thus given by 



where, as usual, in this order of approximation we have neglected the time derivatives and the pressure 
terms. 


The linearization of the metric field equations is a bit more complicated because we need to establish 


suitable gauge conditions. Since the background is Minkowskian, the perturbed Ricci tensor is given by 



( 100 ) 


where hy = — {l/2)6yh^. The term V^Vy4) on the righthand side of Eq. (19) can be combined with 

the terms d^djiy to give the following gauge conditions 



( 101 ) 


With this choice, the linearized equations for the metric become 



( 102 ) 


Eor consistency, to this order Too = p, Tij = 0, T = —p, which leads to 



(103) 


(104) 


Before solving Eqs. (99), (103), and (104), it is worth noting that while the connection equation (14) 


is invariant under constant rescalings of the field 0 —)■ c0, the other field equations do not share this 
invariance. This is manifest in the combinations (1 + 0o) in the above perturbation equations. 


Using the generic solution 



( 105 ) 
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we find that 


hodit,x) = 


Air 


3^ ^—rn^\x—x'\ 




3|a: — x'\ 


K 




47 r(l + 0o) 


K 




d'^x 


X — x'l 

./ p(^, x') 


1 _|_ ^p-m^\x-x'\ 

3 


+ 


Vo |x-X2| 

(1 + 0 o) 6 


X — X' 


-| _ ^ _ Vq |x-X2p 

3 ) (l + <^o) 6 


47 r(l + 0o) 

In spherical symmetry and far from the sources, the above equations become 

K^(t)Q M 

(^(r) = —e 


4o(^) = 


127r r 
2G'efrM , V"o 


r (1 + 0o) 6 ’ 

^ p7G.,Af _ Vi 


where we have denoted 


Geff = 


K 


87r(l + (^o) 


1 + 


^ [1 + 00 exp {-m^r) /3] 

[1 - 0oexp /3]’ 

2 214 — V(j) — 0(1 + 0 ) 10,0 

= -4- 


0=00 


(106) 

(107) 

%.(108) 


(109) 

( 110 ) 
( 111 ) 


( 112 ) 

(113) 

(114) 


These results represent the standard post-Newtonian metric up to second order for this class of theories. 

We emphasize a striking feature of f{X) gravity. Note that in f{R) gravity, to obtain 7 1 from 

(113), there is only one possibility [55,56], namely, m^r S> 1 from millimetres to astronomical scales, 
i.e., the range of the scalar interaction, 1/m^, should be smaller than a few millimetres. In the current 
case, however, there are two possibilities to obtain 7 1 . The first one is the same as in f{R) theories 

and involves a very massive scalar field. The second possibility implies a small value 0o. If 0o 1, 
then the Yukawa-type corrections are very small regardless of the magnitude of m^. This could allow 
for the existence of a long-range scalar field able to modify the cosmological and galactic dynamics, but 
leaving unaffected the Solar System. Subtle modifications could in the most optimistic case be detected 
as anomalies in the local gravitational field [91]. 


4.2. Galactic phenomenology: Stable circular orbits of test particles around galaxies 

The most direct method for studying the gravitational field inside a spiral galaxy is provided by the 
galactic rotation curves. They are obtained by measuring the frequency shifts 2 ; of the 21-cm radiation 
emission from the neutral hydrogen gas clouds. The 21-cm radiation also originates from stars. The 
21 -cm background from the epoch of reionization is a promising cosmological probe: line-of-sight 
velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, 
which linear theory shows can separate cosmological from astrophysical information (for a recent review 
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see [92]). Instead of using z the resulting redshift is presented by astronomers in terms of a veloeity field 
vtg [93,94], 

In the following, we will assume that the gas elouds behave like test partieles, moving in the statie 
and spherieally symmetrie geometry around the galaxy. Without a signifieant loss of generality, we 
assume that the gas elouds move in the galaetie plane 9 = 7r/2, so that their four-veloeity is given by 
= (^, f, 0, 0^, where the overdot stands for derivation with respeet to the affine parameter s. In this 
subseetion we find it illustrative to restore the units of c. 

The statie spherieally symmetrie metrie outside the galaetie baryonie mass distribution is given by the 
following line element 

ds2 = ^ gA ( r)^^2 ^ ^2 ^^ q 2 ^ ^^^^2 ^^^ 2 ^ ^ ^ j ^ 5 ^ 


where the metrie eoeffieients i/(r) and A(r) are funetions of the radial eoordinate r only. The motion of 
a test partiele in the gravitational field with the metrie given by Eq. (115), is deseribed by the Lagrangian 
[95,96] 


L = 



(116) 


where = dO'^ + sin^ Odcf^, whieh simplifies to dVt^ = dcf)^ along the galaetie plane 0 = 7r/2. From 
the Lagrange equations it follows that we have two eonstants of motion, namely, the energy E per unit 
mass, and the angular momentum I per unit mass, given hy E = and I = cr^0, respeetively. 

The normalization eondition for the four-veloeity = — 1 gives 1 = from 

whieh, with the use of the eonstants of motion, we obtain the energy of the partiele as 


^2 


C 


2 


1/+A ’2 \ V 

e r + e 



(117) 


From Eq. (117) it follows that the radial motion of the test partieles is analogous to that of partieles 
in Newtonian meehanies, having a veloeity f, a position dependent effeetive mass meff = 2e^’''^, and an 
energy E. In addition to this, the test partieles move in an effeetive potential provided by the following 
relationship 

+ ( 118 ) 

The eonditions for eireular orbits, namely, dVes/dr = 0 and f = 0 lead to 




^2 

~2 1- rv' 12 ’ 



1 - rz/V2 ’ 


(119) 


respeetively. Note that the spatial three-dimensional veloeity is given by 


t)2(r) 


= e 



( 120 ) 


For a stable eireular orbit dr jdt = 0, and the tangential veloeity of the test partiele ean be expressed as 



= e 


( 121 ) 
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In terms of the eonserved quantities, and along the galaetie plane 0 = 7r/2, the angular veloeity is given 
by 


■'tg 


V) 




and taking into aeeount Eq. (119), we finally obtain the following relationship [95,96] 




rn 

~ 2 ' 


( 122 ) 


(123) 


Therefore, onee the tangential veloeity of test partieles is known, the metrie funetion u{r) outside the 
galaxy ean be obtained as 


i/(r) = 2 


^tg 


[r) dr 


(124) 


The tangential veloeity Vtg/c of gas elouds moving like test partieles around the eenter of a galaxy is 
not direetly measurable, but ean be inferred from the redshift z^o observed at spatial infinity, for whieh 
1 + Zoo = exp [{uoo — n) /2] (1 divtg/c) / — v^g/c^ [95,96]. Due to the non-relativistie veloeities 

of the gas elouds, with vtg/c < (4/3) x 10“^, we observe that vtg/c •zd, Zoo, as the first part of a 
geometrie series. The observations show that at distanees large enough from the galaetie eenter the 
tangential veloeities assume a eonstant value, i.e., Vtg/c ~ eonstant [93,94]. In the latter regions of the 
eonstant tangential veloeities, Eq. (124) ean be readily integrated to provide the following metrie tensor 
eomponent 


e = 


r 

Rv 




l + 2^1n( — 


(125) 


where Ry is an arbitrary eonstant of integration. If we mateh the metrie given by Eq. (125) with the 
Sehwarzsehild metrie on the surfaee of the galaetie baryonie matter distribution, having a radius Rb, 


\r=RB 


= 1 — 2GMb/c^Rb, we obtain the following relationship 

R - 

lL]j 


(1 - ^GMBlc^RBf^"^^ 


(126) 


An important physieal requirement for the eireular orbits of the test partiele around galaxies is that 
they must be stable. Eet tq be the radius of a eireular orbit and eonsider a perturbation of it of the form 
r = tq + 5, where 5 <C tq [97]. Taking expansions of Ks and about r = rg, it follows from 
Eq. (117) that 

5 + (ro) 5 = 0. (127) 

The eondition for stability of the simple eireular orbits requires (rg) > 0 [97]. Henee, with 
the use of the eondition V/g (rg) = 0, we obtain the eondition of the stability of the orbits as 
[SU + ru" > ri^'^/2] \r=ro- By taking into aeeount Eq. (123), it immediately follows stable eireular orbits 
always exist for massive test partieles. 


4.2.1. Galaetie geometry and tangential veloeity eurves in hybrid metrie-Palatini gravity 

The rotation eurves only determine one, namely z/(r), of the two unknown metrie funetions, z/(r) 
and A(r), whieh are required to deseribe the gravitational field of the galaxy [98]. Henee, in order 
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to determine A(r) we proeeed to solve the gravitational field equations for the hybrid metrie-Palatini 
gravitational theory outside the baryonie matter distribution. This allows us to take all the eomponents 
of the ordinary matter stress-energy tensor as being zero. Furthermore, our task beeomes easier when 
we restriet to perturbative weak-field treatment. 

The weak field limit of the gravitational theories at the Solar System level is usually obtained by using 
isotropic coordinates as above in Section 4.1. However, it is useful to apply Schwarzschild coordinates 
in studying exact solutions and in the context of galactic dynamics, and this is what we do here. .We 
assume that the gravitational field inside the halo is weak, so that v{r) ~ A(r) {vtg/cY, which allows 
us to linearise the gravitational field equations retaining only terms linear in {vtg/cY and again consider 
the scalar field as 0 = 0o + where <C 1 is a small perturbation around the background value 0o > 0. 
The Klein-Gordon equation was already solved in Section 4.1 at this limit and the result was that the 
interaction range is given by = 1/m^, where is given in Eq. (114). In a spherically symmetric 
configuration, the general solution then has the form 


ip{r) = ^0- 




(128) 


where tho is an integration constant. Comparing this expression with the results obtained in [28] for the 
weak-field limit (taking into account the transformation from isotropic to Schwarzschild coordinates), 
we find that 

2GMb 




= 


< 0 , 


(129) 


where Mb and are the mass and the radius of the galactic baryonie distribution, respectively. 
Within this linear approximation the stress-energy tensor of the scalar field is given by 




K 


where a and /3 are defined by 


a = \-\v'{4)0) , l3 = -\v{4)o). 


(130) 


(131) 


Therefore the linearized gravitational field equations take the form 

1 / AN 
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- 1 - 2 ~ ^ ’ 

-- I n" + - -^ 1 = aif + {3 . 


Eq. (132) can be immediately integrated to provide 

An . 1 r . . .A An 


/S 


A(r) = — + - {aif + /3) r^dr = — -f — r 




1 + — 


(132) 

(133) 

(134) 

(135) 


where Aq is an integration constant. Comparing again with the results obtained in [28] for the weak-field 
limit, we find that Aq = 2GM/c^. The tangential velocity of the test particles in stable circular orbits 
moving in the galactic halo can be derived immediately from Eq. (133), and is given by 
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(136) 
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which in terms of the solutions found above beeomes 


Vo 2 , GMb 
—- = —r H- 

c 2 g ^ 2 ^ 2 r 


l + f )(2 + arJ) + ^(l + arJ) 


(137) 


where Vq = —/S = V (0o) /2. The term proportional to eorresponds to the eosmologieal baekground, 
namely the de Sitter geometry, and we assume that it has a negligible eontribution on the tangential 
veloeity of the test partieles at the galaetie level. 

On the surfaee of the baryonie matter distribution the tangential veloeity must satisfy the boundary 
eondition 

< (Rb) GMb 


C2 C^B ’ 

whieh, with the use of Eq. (129), gives the following eonstraint on the parameters of the model, 

1 + ^) (2 + + ^(1 + ~ 0 • 

' ^ / 'if 

In order to satisfy the above eondition would require that — 2 < < —1, or, equivalently. 


(138) 


(139) 


E'(0o)>O, 2<-f"'(0o)r2 <3 


(140) 


In the regions near the galaetie baryonie matter distribution, where Rb < r r^, we have e ~ 
1, to a very good approximation. Henee in this region the tangential veloeity ean be approximated as 

v'} 2GMb — c^'^o (ar'^ + 2) + 2 vho . nx 

" ~ ^ ^ ^ ^0^^^-;^(l + «r5)r, i?B<r<r^.(141) 


2c^r 


2r„ 


K 


If the parameters of the model satisfy the eondition 


2GMb - («rj + 2) ^ 0 , 


(142) 


the term proportional to 1/r beeomes negligible, while for small values of ^nd \ar^ 2 1 ~ 1, the term 
proportional to r ean also be negleeted. Therefore for the tangential veloeity of test partieles rotating in 
the galaetie halo we obtain 
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Rb < r r, 
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(143) 


Sinee aeeording to our assumptions, 3> 1, the eoeffieient a ean be approximated as a —V (0o) /2, 
whieh provides for the rotation eurve, in the eonstant veloeity region, the following expression 


< ^oV'{<t2o)G 


Rb < r 


(144) 


Sinee 4/o < 0, the sealar field potential must satisfy the eondition V' (0o) < 0. In the first order of 
approximation, with exp {—r/ry^) ^ 1 — r/r^, for the tangential veloeity we obtain the expression 
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(145) 
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Alternatively, in general we ean write the tangential veloeity as follows. 


-"tg 


V^O 2 , GMb 

—r H- 
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200 

1 + —e 
3 




(146) 


GM-QjeP'—r R^—r 

As eompared to our previous results, in this representation we have e '"‘p instead of e . Sinee 
we are working in a regime in whieh Rb <C r^p, the ehoiee of the eonstants Rb or Mb does not seem very 
relevant, sinee it just amounts to a resealing of 0o. From now on we will also assume that ^ 

From the above equation we want to find the eonstraints on the model parameters that arise from the 
expeeted behavior at different seales. For that purpose, it is eonvenient to write the equation, equivalently, 
as follows: 
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1 + ^(2 + arJ)e 
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(147) 


At intermediate seales, the asymptotie tangential veloeity tends to a eonstant. If we expand the 

_ r 

exponential as e 1 — r/r^, hen we obtain the following three eonstraints on the free parameters of 
the model. 


0 1 + (2 + arl) ^ 0 , 


h) (2 + ar^) ^ C = constant, 
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GM. 
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(148) 

(149) 

(150) 


With inereasing r, and by assuming that the eondition r still holds, the rotation eurves will deeay, 
at very large distanees from the galaetie eenter, to the zero value. 


4.2.2. On astrophysieal tests of hybrid metrie-Palatini gravity at the galaetie level 


In [98], some observational astrophysieal tests of hybrid metrie-Palatini gravity at the galaetie level 
were diseussed. More speeifieally, a eomparison of the theoretieal predietions of the model with a 
sample of rotation eurves of low surfaee brightness galaxies was analysed. Indeed, the predietions of 
the theoretieal model with the observational results show that the eontribution of the sealar field energy 
density to the tangential veloeity of the test partieles ean explain the existenee of a eonstant rotational 
veloeity region around the baryonie matter, without requiring the presenee of the dark matter. The 
possibility of observationally determining the funetional form of the sealar field (p by using the veloeity 
dispersion of stars in galaxies, and the red and blue shifts of gas elouds moving in the galaetie halo 
eould also be worth eonsidering. It ean be shown that one ean eonstrain the explieit funetional form of 
the sealar field, and the free parameters of the model, in order to adequately fit the observational data. 
We refer the reader to [98] for more details. Of eourse, was partiele matter deteeted direetly, we eould 
exelude the gravitational f{X) alternative to dark matter. 
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4.3. Galactic clusters: The generalized virial theorem in hybrid metric-Palatini gravity 

A large number of astronomical and astrophysical observations confirm that galaxies form complex 
hierarchical structures, in which galaxies concentrate in large objects, called clusters of galaxies, 
bounded by the gravitational interaction. The total mass of the galaxy cluster ranges from IO^^Mq for 
the so-called groups (formed by a few hundred of galaxies) up to a few IO^^Mq for very large clusters, 
containing thousands of galaxies. From a morphological point of view galactic clusters are usually by 
a main component, which is regular and centrally peaked [99,100]. For fundamental physics research 
the main importance of the galactic clusters consists in the fact that they are considered to be “dark 
matter” dominated astrophysical objects. Moreover, their formation and evolution is almost entirely 
controlled by the gravitational force, a property which allows the testing of different dark matter models. 
On the other hand the mass distribution inside the clusters is fully determined by the initial conditions 
of the mass distribution that originate in the early universe [101-103], thus allowing the use of galaxy 
cluster properties to also test cosmological models. In this subsection we discuss the possibility of 
observationally testing the metric-Palatini gravity theory by using observational evidence from galaxy 
clusters [104]. 

4.3.1. Galaxy cluster as a system of identical and collisionless point particles 

As a first step necessary to obtain our main result, consisting in the generalization of the relativistic 
virial theorem for galaxy clusters in the hybrid metric-Palatini gravitational we write down the 
gravitational field equations and the Boltzmann for a static and spherically symmetric distribution of 
matter. We adopt a simplified physical model, in which the galactic cluster consists of a self-gravitating 
system of identical, collisionless point particles (the galaxies) in random dynamical motion. The metric 
will thus be described by (115), and the fluid by a distribution function fs that obeys the general 
relativistic Boltzmann equation [105]. 

The energy-momentum tensor of the matter in the cluster, is thus determined by the distribution 
function /s, and its components are given by the equation [106,107] 

= J fBmu^u^ du, (151) 

where m is the mass of the particle (galaxy), = {ut,Ur,UQ,u^) is the corresponding galactic 
four-velocity, with Ut denoting the temporal component. Finally, by du = durduodu^/ut we denote 
the invariant volume element in the velocity space. Alternatively, the energy-momentum tensor 
describing the matter distribution in a cluster of galaxies can be represented in terms of an effective 

ir) 

energy density pes, and of two effective anisotropic thermodynamic pressures, the radial and the 
tangential pressure, respectively. These thermodynamic parameters of the cluster are obtained by 
averaging over the matter and velocity distributions, and are given by 

Peff = = p{ul) = p{ul), (152) 

where p is the mass density of the ordinary baryonic matter, and {uf), i = t,r,9,ip denotes the average 
value of u^, i = t,r, 6, p, representing the square of the components of the four-velocities of the galaxies 
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in the cluster [108]. The full field equations for the metric (115) with the source (152) were listed in Ref. 
[104], but here we shall only need their trace: 




= 47r- 


G 


1 + 
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^ V(4>)+ ^ 
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(2VtV‘ + V„V“) 


—(153) 

0(1 + 0 ) 


where we have denoted = {uf) + {uD + {uD + {ul). 

Since in the following we are interested only in astrophysical applications at the extra-galactic cluster 
scale, we will adopt a Newtonian type approximation, which consists in assuming that the deviations 
from standard general relativity (corresponding to the background value 0 = 0 ) are small for the systems 
we are considering. This approximation implies that 0 -C 1. Thus, Eq. (153) can be approximated as 
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(154) 


where 

AirGp^;^^^ ~ (0) + (2VzV‘ + V„V“) 0 - |Vt0VV, (155) 

corresponds to an effective, geometric type "energy" of the scalar field in the hybrid metric-Palatini 
gravitational model. 


4.3.2. The relativistic Boltzmann equation 


Next we proceed to the second step in the derivation of the virial theorem for galaxy clusters 
and determine the general relativistic Boltzmann equation that governs the evolution of the galactic 
distribution function fs- 

A basic result in statistical physics is the transport equation for the distribution function for a system of 
particles in a curved arbitrary Riemannian space-time. This transport equation is given by the Boltzmann 
equation without collision term, and which can be formulated as [106,107] 

(!>“^-pVr-4)/B = o. (156) 

where p°‘ is the four-momentum of the galaxy (particle), and T^^^ are the Christoffel symbols associated 
to the metric. An important consequence of the collissionless Boltzmann equation is that the local phase 
space density, as measured by an observer in a frame co-moving with a galaxy, is conserved. 

An important simplification of the mathematical formulation of the Boltzmann equation can be 
achieved by introducing an appropriately chosen orthonormal frame, or tetrad el{x), a = 0,1, 2, 3. 
The tetrad fields vary smoothly over some coordinates neighborhood U, and they satisfy the general 
condition g^’''e^e^ = for all x G f/, where 77 “^ denotes the Minkowski metric tensor [106-108]. A 
basic property of the tetrad fields is that any tangent vector defined at an arbitrary point x can be 
represented as p^ = p°‘e^, a relation which defines the tetrad components 

For the spherically symmetric line element given by Eq. (115), the frame of orthonormal vectors can 
be chosen in an appropriate way as [106-108]: 


= el = e^^^5l, el = T5l, el = Tsme6l 


(157) 
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The tetrad eomponents of the four veloeity are = u^e^. In the tetrad eomponents eorresponding to 
our present ehoiees the relativistie Boltzmann equation Eq. (156) is given by 








= 0 , 


(158) 


where the distribution funetion fs = is a funetion of and respeetively, and 76c = 

Rieei rotation eoeffieients [106-108]. Due to the spherieal symmetry of the problem the 
distribution funetion depends only on the radial eoordinate r, and henee Eq. (158) beeomes [108] 
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(159) 


Due to the spherieal symmetry of our astrophysieal system the eoeffieient of cot 9 in Eq. (159) must 
be zero. Erom a mathematieal point of view this implies that the distribution funetion fs is only a 
funetion of r, m and m| + m|. As a next steps in our analysis we multiply Eq. (159) by murdu, and we 
integrate over the veloeity spaee. Then, by taking into aeeount that the distribution funetion fs vanishes 
suffieiently rapidly as the veloeities tend to ±cx), we find the equation 

4 [P (“?>] + \l> [{<) + <“?>] - 2 («?>] = 0. (160) 

Now we multiply Eq. (160) by dvrr^, and by integrating over the entire eluster volume, eorresponding to 
a radius R, gives [108] 

^ 47rp [(m^) + (m^) + (m^)] rMr - \ 47rrV [(“o) + («?>] = 0- 

4.3.3. Geometrieal quantities eharaeterizing galaetie elusters 


In order to obtain some analytieal estimations of the main geometrieal and physieal quantities 
eharaeterizing the galaetie elusters we introduee some approximations to the motion of both test partieles 
in stable eireular orbits around galaxies, as well as to the motion of galaxies in galaetie elusters. As a 
first approximation we assume that u and A are slowly varying funetions of the radial eoordinate r. Then 
in Eq. (154) we ean negleet all the quadratie terms as being negligibly small as eompared to the first 
order terms. Seeondly, astronomieal observations show that the motion of the galaxies in elusters is 
non-relativistie. Henee the galaetie veloeities are mueh smaller than the veloeity of the light, that is, 
(m^) (m|) (uf) <C (mq) ~ 1. Thus, Eqs. (154) and (161) ean be written as 
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is the total kinetic energy of the galaxies in the cluster. We define the total baryonic mass Mb of the 
galactic cluster as 


pR pR 

Mb = / dM(r) = / dvrpr^dr (165) 

Jo Jo 

We assume that the main contribution to the baryonic mass Mb is due to the presence of the intra-cluster 
gas and of the stars. On the other hand we also include in Mb the mass contribution of other particles, 
like, for example, massive neutrinos, which may also give a significant contribution to Mb- 
By multiplying Eq. (162) by r^, and integrating from 0 to r we obtain 

GMB{r) = (r), (166) 

where we have introduced the notation 

(r) = Air ( . (167) 

Jo 

It is interesting to note at this point that in hybrid metric-Palatini gravity, the quantity has 

essentially a geometric origin. Hence it is natural to call it as the geometric mass of the galactic cluster. In 
the following we introduce the gravitational potential energies of the cluster by means of the definitions 
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where R is the cluster radius. By multiplying Eq. (166) with dMB{r), and by integrating from 0 to the 
cluster radius R, we obtain the important relation 
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4.3.4. The generalized virial theorem in hybrid metric-Palatini gravity 


(170) 


As a last step in our analysis, with the help of Eq. (163), we obtain the generalization of the virial 
theorem in hybrid metric-Palatini gravity, which can be formulated in the familiar form 

2K + n = 0. (171) 

In the above equation the total gravitational potential energy of the system, fl, defined as 

n = nB- , (172) 

contains a supplementary term which has a purely geometric origin. 

It is useful to represent the generalized virial theorem, given by Eq. (171), in a more transparent 
physical form, which can be obtained by introducing the radii Ry and R^, defined by 

Rv = Ml/[ ^^M!ldMB(r), (173) 

' Jo 
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and 
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(174) 


respectively. The quantity R^, having a geometric origin similarly to the geometric mass considered 
above, may be called as the geometric radius of the cluster of galaxies in the hybrid metric-Palatini 
gravity theory. Hence, finally we obtain the baryonic potential energy VLb and the effective scalar field 
potential energy as given by 
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(175) 


respectively. Another important observational quantity, the virial mass My of the cluster of galaxies is 
defined as follows 
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GMbMv 

Ry 


(176) 


Eventually, the fundamental relation between the virial and the baryonic mass of the galaxy cluster 
can be obtained after substitution of Eq. (176) into the virial theorem as 
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If My/M b > 3, a condition which is satisfied by the astrophysical parameters of most of the 
observed galactic clusters, then from Eq. (177) we obtain the virial mass of galactic clusters in hybrid 
metric-Palatini gravity, which can be approximated as 
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(178) 


The virial mass My is determined observationally from the study of the velocity dispersion a/: of 
the stars and of the galaxies in the galactic clusters. An important consequence of the virial theorem 
in hybrid metric-Palatini gravity is that in a cluster with mass Mtot most of the mass is in the form of 
the geometric mass M^^^\ Hence we can use the approximation ^ Mtot in order to study the 

cluster dynamics. A fundamental question related to the possibility of the observational testing of the 
hybrid metric-Palatini gravity is to find out what astrophysical processes can detect the presence of the 
geometric mass. Such an observational possibility may be provided by gravitational lensing. Through 
the study of the lensing properties of the galactic clusters one can obtain direct evidence of the existence 
of the geometric mass, of its distribution properties, as well as of the gravitational effects associated to 
the presence of the scalar field. It is interesting to note that gravitational lensing can give us theoretical 
information even at cosmical scales extending far beyond of the virial radius of the matter distribution 
of the galactic cluster. 
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4.3.5. On astrophysical tests of hybrid metric-Palatini gravity at the eluster level 

In eoneluding, in the framework of hybrid f{X) gravity theory we have established the the existenee 
of a striet proportionality between the virial mass of the eluster and its baryonie mass, a relation whieh 
ean also be tested observationally. One of the important, and observationally testable, predietions of the 
hybrid metrie-Palatini gravitational “dark matter” model is that the geometrie masses assoeiated to the 
elusters, as well as its gravitational effeets, extend beyond the virial radii of the elusters. Observationally, 
the virial mass My is obtained from the study of the veloeity dispersions of the stars in the eluster. Due 
to the observational uneertainties, this method eannot give a reliable estimation of the numerieal value of 
the total mass Mb + in the eluster. However, a mueh more powerful method for the determination 

of the total mass distribution in elusters is the gravitational leasing of light, whieh may provide direet 
evidenee for the gravitational effeets at large distanees from the eluster, and for the existenee of the 
geometrie mass. The presenee of hybrid metrie-Palatini modified gravity effeets at large distanees from 
the eluster, and espeeially the large extension of the geometrie mass, may lead to signifieantly different 
leasing observational signatures, as eompared to the standard relativistie/dark matter model ease. The 
bending angle in the hybrid metrie-Palatini gravity models eould be larger than the one predieted by 
the standard dark matter models. Therefore, the detailed observational study of the gravitational leasing 
eould diseriminate between the different theoretieal models introdueed to explain the motion of galaxies 
(“partieles") in the elusters of galaxies, and the standard dark matter models. We refer the reader to Ref. 
[104] for more details. 

Finally, it is worth pointing out that hybrid gravity ean be preeisely tested also at smaller seales like 
those around the Galaetie Centre. As reported in [109], the observed peeuliar orbit of S2 star, moving 
around the eentre of our Galaxy, is theoretieally reeonstrueted if one adopts the gravitational potential 
eoming from hybrid gravity. This result opens new perspeetives in aehieving preeision tests for the 
theory. 

5. Conclusions 

In this work we have presented a hybrid metrie-Palatini framework for theories of gravity, and have 
tested the new theories it entails using a number of theoretieal eonsisteney eheeks and observational 
eonstraints. From the field theory perspeetive, we found that the f{X) elass of theories, where X = 
R -I- k^T, enjoys a similar speeial status amongst the more general hybrid metrie-Palatini theories as the 
f{R) theories within the narrower framework of purely metrie gravity [70]. This is so beeause when 
one exeludes theories inhabited by ghost-like, superluminally propagating and otherwise pathologieal 
degrees of freedom, there is evidenee, as shown in Seetion 2.4, that the f{X) family is singled out as 
the only viable form of an aetion one ean eonstruet using the metrie (and thus the metrie Levi-Civita 
eonneetion) and an independent “Palatini eonneetion”. The underlying reason is that in the speeial ease 
of f{X) aetions the higher derivatives in the gravity seetor ean be separated into a sealar mode, thus 
avoiding an Ostrogradskian instability. Indeed the f{X) gravity represents a generie ease within the 
one-parameter family of the Algebraie Sealar-Tensor theories (reeall Eq. (17)), at one end of whieh lies 
the pure Palatini f{TZ) (wherein the sealar field is a funetion of the stress-energy traee T) and at the other 
end the pure metrie f{R) (where the field is a funetion of the metrie eurvature R). Furthermore, the 
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propagating degrees of freedom have proven to be healthy also on eurved baekgrounds as eonfirmed also 
by our eosmologieal perturbation analysis in Seetion 3.2. Coneeming the Cauehy problem, it was shown 
that in this elass of theories the initial value problem ean always be well-formulated and well-posed 
depending on the adopted matter sourees. 

Having established the theoretieal eonsisteney and interest on the hybrid metrie-Palatini f{X) family 
of theories, we eonsidered applieations in whieh these theories provide gravitational alternatives to dark 
energy. As shown by our post-Newtonian analysis in Seetion 4.1, the hybrid theories are promising 
in this respeet as they ean avoid the loeal gravity eonstraints but modify the eosmologieal dynamies at 
large seales. This is simply beeause as a sealar-tensor theory, the hybrid f{X) gravity is eharaeterised 
by an evolving Brans-Dieke eoupling, whieh allows to introduee potentially large deviations from GR 
in the past (and future) as long as the eoupling at the present epoeh is strong enough to hide the field 
from the loeal gravity experiments. In eontrast, in the metrie f{R) models the Brans-Dieke eoupling 
is a finite eonstant and one needs to invoke some of the various “sereening meehanisms” (workings of 
whieh remain to be studied in the hybrid theories) in order to reeoneile the Solar system experiments 
with eosmology. 

eosmologieal perturbations have been also analysed in these models up to the linear order [28,54, 
78], and the results imply that the formation of large-seale strueture in the aforementioned aeeelerating 
eosmologies is viable though exhibits subtle features that might be deteetable in future experiments. In 
Seetion 3.2, we derived the full perturbations equations and extraeting their Newtonian limit, deseribing 
the observable seales of the matter power speetrum, the growth of matter overdensities was shown to 
be modified by a time-dependent effeetive fifth foree that is expeeted to modify the redshift evolution 
of the growth rate of perturbations. We also note that numerieal studies of the perturbations imply 
that the differenee of the gravitational potentials ean exhibit oseillations at higher redshifts even when 
the baekground expansion and the full lensing potential are indistinguishable from the standard ACDM 
predietions [78]. Sueh features eould potentially be observed in eross-eorrelations of the matter and 
lensing power speetra, but detailed ealeulations of the eosmie mierowave baekground anisotropies and 
other perturbation observables remain to be earried out. This is espeeially worthwhile in view of the 
potential of the fortheoming Euelid mission to experimentally test different aeeelerating eosmologieal 
models by their impaet to the large-seale strueture [110]. 

At an effeetive level, the f{X) modifieations involve both (the traee of) the matter stress energy and 
(the Rieei sealar of) the metrie eurvature, and from this point of view it appears appealing to speeulate on 
the possible relevanee of these theories to both the problems of dark energy and dark matter, in a unified 
theoretieal framework and without distinguishing a priori matter and geometrie sourees. Various aspeets 
of dark matter phenomenology from astronomieal to galaetie and extragalatie seales were diseussed in 
seetion 4. The generalised virial theoreom ean aequire, in addition to the eontribution from the baryonie 
masses, effeetive eontributions of geometrieal origin to the total gravitational potential energy, whieh 
may aeeount for the well-known virial theorem mass diserepaney in elusters of galaxies. In the eontext of 
galaetie rotation eurves, the sealar-field modified relations between the various physieal quantities sueh 
as tangential veloeities of test partieles around galaxies, Doppler frequeney shifts and stellar dispersion 
veloeities were derived. More reeently, observational data of stellar motion near the Galaetie eentre was 
eompared with simulations of the hybrid gravity theory, whieh turned out partieularly suitable to model 
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star dynamics. Yet, to promote the /(X) theory into a eonvineing alternative to partiele dark matter, one 
should produee also the eosmologieal sueeesses of the ACDM model without the CDM eomponent. 

Though wormhole solutions have already been found in these theories [111], the nature of possible 
blaek hole solutions remains an outstanding open question. Though no-hair theorems for sealar-tensor 
theories exelude the simplest nontrivial solutions, their assumptions are restrietive. Another interesting 
though yet unstudied issue is the strong field regime of hybrid gravity and the eonstraints that one ean 
put on the theories from astrophysieal data sueh as measurements of binary pulsars. Finally, the possible 
relevanee of the hybrid gravity framework to the early universe eosmology has not been explored in 
any detail. The Einstein frame formulation of the sealar-tensor deseription (25) eould provide a useful 
starting point to investigate how the inflaton potential ehanges due to finite VLa, in order to understand 
how the ’’hybrid” nature of spaeetime modifies the standard predietions of e.g. the 'R? model of inflation. 

To eonelude, whilst the physies of the metrie and the Palatini versions of f{R) gravity have been 
uneovered in exquisite detail in a great variety of different eontexts [9-18], those studies largely wait to 
be extended for the hybrid f{X) version of the theory. We believe the results this far, as reported in this 
review, provide eompelling motivation for the further exploration of these partieular theories. 
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